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Abstract. The possibilities of observing some nonlinear electrodynamic effects, which can be manifested in 
hard emission of X-ray, gamma ray pulsars and magnetars by X-ray and gamma ray astronomy methods are 
discussed. The angular resolution and sensitivity of modern space observatories give the opportunity to study 
the nonlinear electrodynamic effects, which can occur in very strong magnetic fields of pulsars (B ~ 10 12 G) and 
magnetars (B ~ 10 G). Such magnetic field magnitudes are comparable with the typical value of magnetic field 
induction necessary for manifestation of electrodynamics non-linearity in vacuum. Thus, near a magneticneutron 
star the electromagnetic emission should undergo nonlinear electrodynamic effects in strong magnetic fields (such 
as bending of rays, fluxes dispersing, changing of spectra and polarization states). Manifestations of these effects 
in detected hard emission from magnetic neutron stars are discussed on the base of nonlinear generalizations 
of the Maxwell equation in vacuum. The dispersion equations for electromagnetic waves propagating in the 
magnetic dipole field were obtained in the framework of these theories. The possibility of observing the bending 
of a ray and gamma ray flux dispersing in the neutron star magnetic field are analyzed. The only nonlinear 
electrodynamicseffect, which can be measured principally, is the effect of gamma ray flux dispersion by the 
neutron star magnetic field. Studying this effect we can also obtain information on the nonlinear electrodynamics 
bending of a ray in the source. The main qualitative difference in predictions of different nonlinear electrodynamics 
theories are discussed. 



1. Introduction 

The advanced X-ray and gamma-ray space observatories, 
such as the well known Chandra (Van Spreybroeck et al. 
1997; CXC 2000 x ), XMM-Newton (HEASARC 2001 2 ), 
and INTEGRAL (Gehrels et al. 1997; Winkler 1999) per- 
mit to observe astrophysical sources in their hard emission 
at the sensitivity levels ~ 10 -5 phot sm _2 s _1 keV -1 (for 
0.1-10 keV photons) and ~ 10~ 6 phot cm^s" 1 keV -1 
(for 0.1-1.0 MeV photons) for real time exposures (~ 10 6 
s) with the angular resolution of ~ 0.1 arc sec (soft X-rays) 
and about ten arc min (gamma-rays). The space missions, 
which are in formulation now, such as Constellation X, 
XEUS, NeXT will improve the sensitivity by 15-100 times 
and MAXIM Pathfinder gives a 1000 times increase in the 
angular resolution (Weaver, White, & Tananbaum 2000; 



1 Chandra MSFC, and Chandra IPI Teams 2000, Chandra 
Proposers Observatory Guide (POG) Electronic Version 2.0 is 
available at: http://asc.harvard.edu/udocs/docs/docs.html 

2 High Energy Astrophysics Science Archive Research 
Center 2001, "A Portable Mission Count Rate 
Simulator (W3PIMMS)", Version 3.1c is available at: 
http:/ /heasarc.gsfc.nasa. gov. /TOOLS/w3pimms. html 



White & Tananbaum 2000; Parmar et al. 1999; Cash, 
White, & Joy 2000). This gives us an opportunity to study 
the physical conditions in very strong gravitational and 
electromagnetic fields near compact relativistic objects. 

In particular, the studies of phenomena in the vicinity 
of a neutron star make it possible to obtain information on 
the properties of matter in the states, which are unattain- 
able in the ground laboratories. In this sense, a problem 
of a great importance is to search and study the vacuum 
nonlinear electrodynamics effects, which can occur in very 
strong magnetic fields of such objects as pulsars (B ~ 10 12 
G) and magnetars [B ~ 10 15 G). 

As it is well known the nonlinear electrodynamics of 
physical vacuum during a long time has no experimental 
confirmation, and was usually regarded as an abstract the- 
oretical model. But now its status has changed drastically. 
The last experiments on light-to-light scattering made at 
Stanford (Burke et al., 1997) show that electrodynamics 
in vacuum is really a nonlinear theory. Thus, the different 
models of nonlinear electrodynamics of the vacuum and 
their predictions (Ginzburg 1987; Alexandrov et al. 1989; 
Rosanov 1993; Bakalov, D. et al. 1998, Denisov 2000a,b; 
Rikken & Rizzo 2000, Denisov & Denisova 2001b, c) should 
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be seriously tested in experiments. However, the magnetic 
fields (B ~ 10 6 G) available in ground laboratories give no 
such opportunity because theory predicts that the typical 
value of magnetic field induction necessary for essential 
manifestation of electrodynamics nonlinearity in vacuum 
is B ~ 4.41 • 10 13 G. 

Since magnetic fields of some pulsars can be charac- 
terized by such magnitudes, and for magnetars can reach 
much greater values, this permits us to conclude that 
nonlinear effects of electrodynamics in vacuum should be 
most pronounced in the vicinity of such astrophysical ob- 
jects. Near a magnetic neutron star electromagnetic emis- 
sion undergoes the nonlinear electrodynamics effects of 
strong magnetic fields. As a result, electromagnetic rays 
are bended, the emission fluxes are dispersed and their 
spectra and polarization states change. The presence of 
super strong magnetic field furthers the forming of a quite 
extended magnetosphere with the radius of about several 
radii of a neutron star. This magnetosphere is opaque, as 
a rule, for the low-frequency part of electromagnetic spec- 
trum. It will be transparent only for X-rays and gamma 
rays. Thus the spectrum, polarization and other parame- 
ters of the detected hard emission from magnetic neutron 
stars are the only data, which can give us the information 
on the main regularities of the nonlinear electrodynamic 
interaction of such emission with the strong magnetic field 
of a pulsar or magnetar. 

Thus, the most appropriate astrophysical objects 
where the nonlinear electrodynamics effects can be man- 
ifested more clearly are some types of rotation-powered 
pulsars, accretion-powered pulsars and magnetars. 

Different nonlinear electrodynamic effects in the vicin- 
ity of a strongly magnetized neutron star were previously 
studied in the context of quantum electrodynamics. In 
particular, vacuum birefringence and its effect on the emit- 
ted spectra and on the propagation of photons in the 
neutron star magnetosphere was discussed in the book of 
Meszaros (1992). It was found that vacuum effects dom- 
inate the polarization properties of the normal modes of 
the near-neutron star medium. This gives rise to a signifi- 
cant change in the medium opacity, thus the polarization 
properties and transport of X-ray radiation from a neu- 
tron star's magnetosphere can be altered by the magnetic 
vacuum effects (Meszaros & Ventura 1978; Meszaros & 
Ventura 1979; Borner & Meszaros 1979; Meszaros et al. 
1980; Meszaros & Bonazzola 1981, Denisov et al. 2002). It 
was also shown that magnetic vacuum effects can change 
the spectra of emitted radiation which leads to a signature 
in the spectra of x-ray pulsars (Ventura et al. 1979). The 
nonlinear quantum electrodynamic effects induced by an 
nonhomogeneous and non-stationary magnetic field of a 
neutron star, including the light bending in the plane of 
the magnetic dipole equator, photon pair production and 
the frequency doubling and modulation at the scattering 
of low frequency electromagnetic waves by the magnetic 
field of an inclined rotator were discussed by Gal'tsov and 
Nikitina, 1983. The quantum electrodynamic effects in the 
accreting neutron stars, in particular, one and two-photon 



Compton scattering in strong magnetic field and its effect 
on the radiation processes (Bussard et al. 1986) as well as 
the vacuum polarization effects in the field of a charged 
compact object (De Lorenci et al., 2001) were also stud- 
ied. However, the analysis presented above concentrated 
on the validity of quantum electrodynamics without com- 
parison with possible alternative theories. Thus, we try to 
obtain some specific predictions by using post-Maxwellian 
items of different nonlinear generalizations of electrody- 
namics and then compare the predictions of different the- 
ories with the goal to use the nonlinear electrodynamics 
effects in neutron stars as a test for post- Maxwellian ef- 
fects. 

2. Nonlinear models of vacuum electrodynamics 

As it is well known, Maxwell electrodynamics is a linear 
theory in the absence of matter. Its predictions concern- 
ing a very wide field of problems (except the subatomic 
level) are constantly confirming with better and better 
accuracy. Quantum electrodynamics, which is based on 
Maxwell electrodynamics complemented by the rcnormal- 
ization procedure, also describes with good accuracy the 
various subatomic processes, and, according to common 
opinion, is one of the best physical theories. 

However, some fundamental physical reasons indicate 
that Maxwell electrodynamics is only the first approxima- 
tion of more general nonlinear vacuum electrodynamics, 
which can be applied in the limit of weak electromagnetic 
fields. 

The electromagnetic field equations, which can be ob- 
tained using the Lagrange formalism, in any nonlinear 
model of vacuum electrodynamics are equal to: 



1 <9D 

curl H = -— — , div D = 0, 

c at 

curl E = — , 

c dt 



(1) 



div B = 0. 



However, in these equations the vectors D and H de- 
pend on vectors B and E, differently in various models, 
since they are defined by the different dependencies for 
the Lagrangian L = L(B, E): 



D = 4tt 



dL 

8E' 



H = -4tt 



dL 



(2) 



At the present time several nonlinear generalizations 
of the Maxwell equations in vacuum are considered in the 
framework of the field theory. The most well known among 
them are the Born-Infeld (BI) electrodynamics (Born & 
Infeld 1934) and the Heisenberg-Euler (HE) electrody- 
namics (Heisenberg & Euler 1936). These theories are 
based on absolutely different principles, and, as a result 
they lead to different electromagnetic field equations. 

2.1. The Born-Infeld nonlinear electrodynamics 

Born and Infeld in their research proceeded from the idea 
of a limited value of the electromagnetic field energy of a 
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point particle. This and some other reasons led them to 
the following Lagrangian of the nonlinear electrodynamics 
in vacuum: 



Ana 2 



y/l + a 2 (B 2 - E 2 ) - a 4 (B • E) 2 - 1 



(3) 



where a - is the constant with units reciprocal to the units 
of magnetic field induction. 

Born and Infold defined this constant from the assump- 
tion that the origin of all rest energy of an electron £ is 
electromagnetic. As a result, the estimate 1/a — 9.18T0 15 
G was obtained (Born & Infeld 1934) from the atomic 
physics constraints. 

Thus, although this theory has a definite Lagrangian, 
it is to a large extent phenomenological, and to verify it is 
necessary, first of all, to measure experimentally the value 
of a 2 parameter or at least to estimate its upper limit. 

In view of relations (2) and (3), in the BI nonlinear 
electrodynamics vectors D and H, are the following func- 
tions of vectors B and E : 



D 



H 



E + a 2 (B • E)B 



v / l + a 2 (B 2 -E 2 )-a 4 (B-E) 2 ' 



B - a 2 (B • E)E 



y/l + a 2 (B 2 - E 2 ) - a 4 (B • E) 2 



For the field magnitudes, which can be achieved in the 
ground laboratories, the values a 2 E 2 and a 2 B 2 are much 
less than one. In this case Lagrangian of the BI nonlinear 
electrodynamics can be expanded into the small parame- 
ters a 2 E 2 << 1 and a 2 B 2 << 1 : 



- — (B 2 

8tt 



E 2 ) + — (B 2 

' 32tt L V 



E 2 ) 2 + 4(B-E) 2 . (4) 



The first term of this expansion is the Lagrangian of the 
Maxwell electrodynamics, and the other term is the non 
linear correction to it, which is proportional to the above 
mentioned small parameters. 

It was shown by Cecotti & Ferrara (1987), Denisov 
(2000a), Denisov et al. (2000), that the BI electrodynamics 
has a number of very interesting properties and in many 
ways it is remarkable theory. 

First, as it was mentioned above, the energy of the 
electromagnetic field of a point charge is a finite quantity 
in the framework of this theory. 

Second, the ideology of this theory is very close to 
Einstein's idea of introducing a non-symmetric metric ten- 
sor dk / Gki with the symmetric part corresponding to 
the usual metric tensor and the antisymmetric part, 
corresponding to the electromagnetic field tensor Fik ■ 
Gik = gik + aFik- Using the relations of tensor algebra 
(Denisova and Mehta 1997) it is not difficult to show that 
the Lagrangian (3) can be written as 



47ra 2 

Besides, though the velocity of an electromagnetic 
wave depends on the values of the fields B 2 and E 2 in this 



theory, it does not exceed the speed of light c in Maxwell's 
electrodynamics. It should be noted also that BI electro- 
dynamics can be obtained from more general sypersym- 
metric theories. 

Thus, the BI electrodynamics in many respects consti- 
tutes a distinguished theory. However, to the present time 
this theory is not developed enough because of the lack 
of quantitive estimations of different effects. In particular, 
none of the scientific publications contain any calculations 
in the BI theory, which could give estimates of the prob- 
ability of e~e+ pair production in the SLAC experiments 
(Burke et al., 1997). On the other hand, it is necessary to 
note, that to the present time there are no experiments, 
which would rejected this theory. 

2.2. The Heisenberg-Euler nonlinear electrodynamics 

As it is well known the HE nonlinear electrodynamics is 
based on the quantum electrodynamic (QED) effect of 
electron-positron vacuum polarization by electromagnetic 
fields. 

Thus, the linear Maxwell electrodynamics is only the 
first approximation of a more general nonlinear electro- 
dynamics (in vacuum), which can be used in the case of 
weak electromagnetic fields, when its magnitudes B and 
E are much smaller than the characteristic quantum elec- 
trodynamic value B q = m 2 e? /eh = 4.41 • 10 13 G, where m 
is the mass of an electron, e is the module of its charge, h 
- is the Plank constant. 

The accurate form of the Lagrangian in this theory has 
not been defined yet. However, for "weak" electromagnetic 
fields corrections to the Maxwell Lagrangian in the first 
non- vanishing order of the quantum electrodynamics per- 
turbation theory have a strictly defined form. As it can be 
seen from calculations (Heisenberg & Euler 1936), if elec- 
tromagnetic fields are not strong (B « B q , E « B q ) 
the first two terms in the vacuum electromagnetic field 
nonlinear Lagrangian expansion in the small parameters 
(B 2 - E 2 )/B 2 and (B • E)/S 2 should have the form: 



8tt' 



[B 2 -E 2 



360tt 2 B 2 



|(B 2 -E 2 ) 2 +7(B-E) 2 },(5) 



where a = e 2 /he w 1/137 is the fine structure constant. 

The vectors D and H, which are contained in equations 
(1), are also nonlinear functions of the vectors B and E 
in this theory. In the first non-vanishing approximation of 
quantum electrodynamics their form is: 



D = E + 



H = B + 



a 



457TB 2 



a 



45ttB 2 



|2(E 2 - B 2 )E + 7(B • E)B j, 
|2(E 2 - B 2 )B - 7(B • E)E}. 



Comparing expressions (4) and (5), it is easy to see that 
they can not be reduced to each other by any choice of the 
a 2 constant. This means, that nonlinear electrodynamics 
with Lagrangians (3) and (5) are essentially different theo- 
ries. Thus, experimental verification of predictions of these 
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theories as well as the solving of problem of their adequacy 
to reality are of great interest. 

2.3. Other models of nonlinear electrodynamics 

Other models of nonlinear vacuum electrodynamics are 
also discussed in the field theory. It is quite natural, that 
in other theoretical models of nonlinear electrodynamics 
the coefficients at the terms (B 2 - E 2 ) 2 and (B • E) 2 in 
the Lagrangian expansion can be absolutely arbitrary. 

Thus, to choose that nonlinear electrodynamics, which 
is most adequate to nature, it is necessary to calculate the 
nonlinear effects in different theories and to compare their 
predictions with the results of the corresponding experi- 
ments. 

To make such calculations easier in the approximation 
of a weak electromagnetic field, we will use a parame- 
terized post-Maxwell formalism, which was elaborated by 
Denisov & Denisova (2001a,c). This formalism is similar, 
in some sense, to the parameterized post-Newton formal- 
ism in the theory of gravitation (Will 1981), which is com- 
monly used for calculating different gravitational effects in 
the weak field of the Solar system. 

We will assume, that the main prerequisite for this 
formalism is that the Lagrangian of nonlinear electrody- 
namics in vacuum is an analytical function of invariants 
Ji = (E 2 - B 2 )/B 2 and J 2 = (E • B) 2 /^, at least, near 
their zero values. Thus, in the case of a weak electromag- 
netic field (Ji << 1, J2 << 1) this Lagrangian can be 
expanded into a converging set in integer powers of these 
invariants: 

i^EEw (6) 

n— m— 

Since, at J\ — > 0, J2 — > the theory with Lagrangian 
(6) should be reduced to Maxwell electrodynamics, then 
ioo = 0, L10 = 1. 

For such an approach, a quite definite number of post- 
Maxwell parameters L nm will correspond to each nonlin- 
ear electrodynamics. From the point of view of the exper- 
iments in a weak electromagnetic field, we can conclude, 
that one nonlinear electrodynamics will differ from the 
other only by the values of these parameters. 

If we limit oneself only to a few first terms in expansion 
(6), then according to the parameterized post-Maxwell 
formalism the generalized Lagrangian of the nonlinear vac- 
uum electrodynamics in the case of weak fields can be 
represented as (Denisov & Denisova 2001a, c): 

L = ^{[E 2 - B 2 ] + £[?7i(E 2 - B 2 ) 2 + 4r? 2 (B • E) 2 ]}, (7) 

where £ = 1/-B 2 , and the value of the dimension- less post- 
Maxwell parameters 771 and 772 depend on the choice of the 
model of nonlinear vacuum electrodynamics. 

In particular, in the nonlinear HE electrodynamics 
parameters 771 and 772 have quite definite values 771 = 
a/(457r) = 5.1 • 10~ 5 , 7] 2 = 7a/(1807r) = 9.0 • 10~ 5 , while 



in the BI theory they can be expressed through the same 
unknown constant a 2 : 771 = i]2 = a 2 L? 2 /4. 

Substituting Lagrangian (7) into expressions (2), we 
obtain vectors D and H of the parameterized nonlinear 
vacuum electrodynamics: 

D = E + 2£{?7i(E 2 - B 2 )E + 2t7 2 (B • E)b}, 

H = B + 2c[t7i(E 2 - B 2 )B - 2t7 2 (B • E)E j. 

Thus, the post-Maxwell formalism without focusing on the 
details of one or another nonlinear electrodynamics, its 
equations, hypotheses and postulates, on all its theoret- 
ical composition takes into account only the final result: 
expansion of Lagrangian, which according to the given 
theory is valid in the weak electromagnetic field approach. 
Further analysis of the theories and revealing of the con- 
cordance of their predictions with experimental results is 
quite general and can be reduced to obtaining the answers 
to two questions: what are the values of post-Maxwell pa- 
rameters in the studied theory and what are the parameter 
values according to the results of corresponding experi- 
ments. 

Thus, one of the goals of this formalism is the calcula- 
tion of "weak" nonlinear electrodynamic effects disrespec- 
tively to any nonlinear theory. The goals of the theory and 
the experiment in this case should be not only the search 
for such an effect (which can refute one or another nonlin- 
ear electrodynamics), but also experiments with the pur- 
poses of determination (with necessary accuracy) of the 
all post-Maxwell parameter values. 

3. The effect of nonlinear-electrodynamic 
bending of a ray 

The electromagnetic emission is the main channel carry- 
ing information on nonlinear electrodynamic effects, which 
can occur in the magnetic dipole field of astrophysical ob- 
jects. The electromagnetic ray is exactly the agent, which 
passing through the neutron star magnetic field under- 
goes nonlinear electromagnetic influence from this field 
independently of the spectral range. Studying the main 
parameters of incoming electromagnetic emission, such as 
dependence of a ray bending angle on impact distance, 
the law of emission intensity decreasing in the course of 
time, etc., it is possible (Denisov, V.I. et al. 2001) to re- 
veal the main dependencies of nonlinear electrodynamic 
interactions of electromagnetic fields. 

3.1. Dispersion Equation 

In order to study the laws of weak electromagnetic waves 
propagation in the dipole magnetic field of a neutron 
star, we will obtain the dispersion equation. We will as- 
sume that a "weak" plane electromagnetic wave propa- 
gates through the permanent magnetic field B of a neu- 
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tron sta. Then in the geometric optics approach we can 
write the following relations: 



E = e cxp[— i(ojt — k • r)], 



(8) 



B = B 



b exp[— i(u>t — k • r)], 



where uj is the frequency, k is the wave vector and vec- 
tors b and e are slowly changing functions of t and r, in 
comparison with the exp[i(uit — k • r)] function. 

Under this approach the dispersion equation can be 
obtained from Lagrangian (7) directly. However, though 
the obtained result will be true, quite legitimate questions 
about its correctness will arise during calculations. Thus, 
to ensure the necessary accuracy of calculations, we will 
add to Lagrangian (7) the terms of higher approximations 
and, hence, write it with surplus accuracy: 



8tt{ 



E 



B 2 ] + S[ m (& - B 2 ) 2 + 4t7 2 (B • E) 2 ] + (9) 



+e 2 [%(E 2 - B 2 ) 3 + r, 4 (E 2 - B 2 )(B • E) 2 ]}. 

This relation is used to obtain the dispersion equation. 
The details of the transitions are presented in Appendix 
A. The result is that according to equations of the post- 
Maxwell nonlinear vacuum electrodynamics in the pres- 
ence of a permanent and regular magnetic field " weak" , 
generally, plane electromagnetic waves of two types can 
propagate in any direction. The corresponding dispersion 
equations are: 



wi(k) - ck[l -M[ kx Bo] 2 + 0(£ 2 Bg)}, 



(10) 



W2(k) = cfc{l -^[kx B ] 2 + 0(£ 2 B*)}. 



It is necessary to note, that the same dispersion equations 
can also be obtained from a simpler Lagrangian (7). 

As it was shown previously (Denisov 2000a) , the exact 
dispersion equation for electromagnetic wave propagating 
in the magnetic field Bq in the BI theory has the form 



(l + a 2 B 2 )^--k 2 -a 2 (k B ) 2 = 0, 



(11) 



independently of its polarization at the any a 2 B 2 value. 

The solution of Maxwell equations (1) for electromag- 
netic waves propagating in the magnetic field shows that 
at r]i 7^ r]2 the waves of both types with the dispersion 
equations (10) are polarized linearly in mutually normal 
planes and propagate with different group velocities. This 
property of electromagnetic waves is well known as bire- 
fringence. 

At r)i — r\2 the both types of electromagnetic waves 
will coincide to the accuracy of terms proportional to £ 2 . 
As a result, electromagnetic waves of the same type with 
arbitrary polarization will propagate in each direction. 

Let us now find now the cikonal equation for an elec- 
tromagnetic wave propagating in the dipole magnetic field 
of a neutron star under the laws of nonlinear vacuum elec- 
trodynamics. For this purpose we will raise relations (10) 



to the second power. Retaining terms linear in £ and tak- 
ing into account that uj = dS/dt, k = VS, we obtain: 

2 / \ 2 



-47/^B 2 ] (Vs) -^(Bo-VS) -0,(12) 

1 - 477 2 ^B 2 ] (Vs)' - 4r7 2 £(Bo-Vs) 2 = 0. 

In the BI theory, as it follows from the relation (11), the 
eikonal equation valid for any a 2 Bg, values has the form 



c 2 V 3 1) 



(13) 



1(^) 2 [1 + a 2 B 2 ] - (VS) 2 - a 2 (VSB ) 2 = 0. 

The solution of these equations in common case is not 
known. Thus, further on we will consider solution of the 
equations (12) only for the rays laying in the dipole mag- 
netic field equator plane. Determination of the evident 
form of eikonal equations (12) for electromagnetic wave 
propagating in the field of a magnetic dipole permits us 
to study the main nonlinear electrodynamic effects, which 
should occur in the magnetic fields of pulsars and magne- 
tars. 

It is necessary to note, that the effects discussed be- 
low can partially be caused by the gravitational field of 
neutron stars. However, in the first order on perturbation 
the gravitational and nonlinear electrodynamic parts are 
additive. Since the gravitational effects were studied re- 
peatedly (Epstein & Shapiro 1980; Will 1981; Meszaros & 
Riffcrt 1988; Riffert & Meszaros 1988) and to the present 
time studies of gravitaional lensing effects in astrophysics 
remains the main goal of some observational programs 
(Sutherland et al. 1996), we will not consider them here 
and will pay attention mainly to the effects of nonlinear 
vacuum electrodynamics. 



3.2. The Bending of a ray from a source located at a 
limited distance from a neutron star in its 
magnetic field 

Let us denote the plane normal to the magnetic dipole 
momentum vector m, as plane XOY. In this case only 
one component of the vector m = |m|e z will be nonzero 
and vector B in this plane can be represented as: B = 
-|m|e z /r 3 . 

Hence the first of the eikonal equation (12) for electro- 
magnetic wave polarized in the XOY plane, which ray lay 
in the same plane, will be: 

1 /dS!\ 

^Vdt 



1- 



4r;i£m 2 



LV dr ) r 2 \ dip 



Oip ) 



0.(14) 



A similar equation with the replacement of the 771 param- 
eter by the 772 parameter and Si by S2 , can be written for 
a ray of electromagnetic wave polarized along the z axis. 

As it is accepted in theoretical mechanics (Landau & 
Lifshitz 1984), we will find the partial solution of equation 
(14) using the variables separation method. As a result, 
we obtain: 



Si = -£ Q t + cup + 



1 



4?7i£m 2 



( 15 ) 
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where £q , a are the constants of integration and all calcu- 
lations were made with accuracy, linear in the small value 
?7 1 ^m 2 /r 6 . 

It should be noted, that in the magnetic equator plane 
the expression (15) is also the solution of BI exact eikonal 
equation (13), if we take into account, that in this theory 
Vi=V2= a 2 B 2 q /4. 

Using relation (15) we can determine the kinematic 
and dynamic parameters describing photon propagation 
in the dipole magnetic field. 

Let us consider the case, when the gamma ray source is 
located at a limited distance l\ from a neutron star or even 
in its nearest vicinity. This can take place in the case of 
an accretion-powered pulsar. The latter conditions can be 
realized in the case of a rotation-powered gamma- pulsar 
(if the polar cap models are valid) . 

Let us denote the distance from neutron star to the 
detector as l^. Then the distance l\ is much smaller than 
I2 and comparable with the neutron star radius R. Hence, 
the dependence of impact distance on time for a circular 
orbit in the first approximation can be represented as: 
b(t) = bo + Ri cos fit, where R\ is the orbit radius, is 
the orbital frequency. 

If we are considering the propagation of a X-ray or 
gamma ray photon from a source located near a galactic 
neutron star, it is convenient to direct the X and Y axes in 
such a way, that a ray from the source travels along the X 
axis with the impact distance b, the center of the dipole 
magnetic field is placed in the center of the coordinate 
system (see Fig. 1) and the spacecraft with the detectors is 
located at the distance I2 from the center of the coordinate 
system near the point x w l 2 . 

Since the value £m 2 /r 6 is small, in order to find the 
bending angles Sipl 2 we can use the algorithm (Darwin 
1961), well established for calculations of angles of light 
gravitational bending. The details of calculations of the 
bending angle 8<p\ in the case of an electromagnetic wave 
polarized in the XOY plane, and 5ip2 for an electromag- 
netic wave polarized along the z axis, are presented in 
Appendix B. For the case, when b < h « h, and l 2 » b 
we obtained following relations for the bending angles: 



Sipl = arctan 



bQi 



+ 



157r?7i£m 2 
466 



(16) 



r , 1577i£m 2 -, . 1 , 



■{ sm[4arcsm( 



166 6 



5ip* 2 



giVJ -16sin[2arcsin(— )]}, 



arctan 



bQ 2 



157r?72^m 2 



r , 15ri 2 £m 2 -, . 1 , 

772^ m 2 r 1 1 , 

-^^{sin[4arcsin(— )] - 16sin[2arcsin(— )]}, 



where 



Qi 



b [ b e - 



_ 1 r 2?7 2 gm 2 



The plus sign in this relation shows, that the magnetic 
dipole field in the magnetic equator plane effects the elec- 
tromagnetic waves as a convex lens. To illustrate depen- 
deces (16) in BI and HE theories we plot 8<p\ and 6ip 2 
versus bjl\ for different values Bo/B q (Fig. 2 - Fig. 4). 

Thus, nonlinear models of vacuum electrodynamics 
with 771 7^ 772, predict different angles of ray bending for 
electromagnetic waves with different polarization. 

It should be noted, that besides the nonlinear elec- 
trodynamic bending (16) of electromagnetic rays will also 
undergo the well known gravitational bending. However, 
because of the different bending angle dependence on the 
impact distance b (1/6 and 1/6 2 in the case of gravitational 
bending (Epstein & Shapiro 1980; Will 1981; Meszaros & 
Riffert 1988; Riffert & Meszaros 1988) and l/b 6 in the 
case of nonlinear electrodynamic bending) mathematical 
processing allows to resolve each of these parts from the 
observational data if the time dependence of impact dis- 
tance b is harmonical b(t) = b + Ri cos fit. 

4. The effect of gamma-ray flux scattering in the 
neutron star magnetic field 

The nonlinear electrodynamic bending of rays in the 
magnetic dipole field differs significantly from gravitation 
bending of these rays. As it is well-known, the ray bending 
in the Schwarzschild gravitation field occurs in the plane 
containing this ray and the star center. Thus, gravitation 
ray bending leads to the increase of electromagnetic flux 
detected by the instrument located on the other side of 
gravitation center. 

In the general case the nonlinear electrodynamic bend- 
ing of electromagnetic ray is not planar. It can be discribed 
as the bending on the series of magnetic force lines, when 
a ray undergoes bending on every force line in the plane 
containing the normal to the surface [Bo k] 2 /fc 2 = const. 
As a result the nonlinear interaction of electromagnetic 
wave with the dipole magnetic field leads in the general 
case to the appearing of the ray curvature and twisting. 

For a complete and detailed study of the principles of 
nonlinear electrodynamic lensing, it is necessary to solve 
Eqs. (12) in the very general case, but not only for rays lay- 
ing in the neutron star magnetic equator plane. However, 
to the present time the solution of these equations is un- 
known. Formulas (16) allow us to study only one particu- 
lar case of nonlinear electrodynamic lensing, when the ray 
source is near the neutron star, at the distance l\ ~ R, 
and electromagnetic emission is a bunch of rays laying in 
the small vicinity of magnetic equator plane. 

When studying gravitation lensing of light emission, it 
can usually be assumed that the distance between light 
source and the gravitating star l\ is larger than the star 
radius R. In this case ray bending leads to the abrupt 
cutting of the field of shadow behind the opaque star (see 
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Fig. 1. The ray bending and shadow region in the case h » R (a); the ray bending and shadow region in the case 
h>R (b). 



• » I f « W » 



F 



■ i ■ i ■ ■ 



/ 



Fig. 2. The BI dependence Sip versus h/b for different values & 2 {b)/B 2 q : (a) B 2 (b)/B 2 = 10~ 2 , (6) B 2 (6)/B 
10, (c) B 2 (&)/£? 2 = 1 (B 2 (fe) is the square of magnetic field at the distance b from the neutron star, 5ip\ = Sip^)- 



Fig. la). As a result, the light flux detected at the large 
distance will be much larger than in the case of gravitating 
center absence (the effect of gravitational lensing) . 

The sources of gamma-rays are located mainly not far, 
or according to some models (Harding, 2000), even near 
the pulsar's or magnctar's surface: l\ <~ R. Thus, even 
with account for the gravitational and nonlinear electro- 
magnetic bending of gamma-rays, we should assume the 
existence of some shadow cone behind the neutron star, 
where gamma-rays never penetrate (see Fig. lb). 

The main condition of the existence of such cone is the 
non-equality: S(p max < arcsin(i?/Zi). Let us consider this 
case in details. 

We will assume, that a magnetized neutron star is lo- 
cated in the coordinate center, the gamma ray source is 
placed in the XOY plane at the distance of l\ from the 
neutron star (left panel of Fig. lb) and the spacecraft car- 
rying the gamma ray detector is located at a large distance 
from the star, at the point xu = h >> h, Vd = p- Let 
us calculate the value of gamma-ray flux I out detected by 



the presence of gravitation and nonlinear electrodynamic 
bending of rays. 

Let us assume, that in the case of the absence of such 
bending, the detector will measure the energy flux Iq = 
A/a[p 2 + (l\ + h) 2 ], where a is the value of gamma-ray 
beam solid angle, A is some constant characterizing the 
gamma-ray source intensity. 

It is quite appropriate to chose the sector of a ring 
of radius, 6, width, db and angle opening 89 « 1 as 
the input aperture. This aperture area will be equal to: 
S in = bdbSe. 

In the input aperture plane the gamma-ray flux 7j n will 
be: I in = A/a[l 2 -R 2 } = I Q [p 2 + (h+l 2 ) 2 ]/[l 2 -R 2 }. Since in 
the field of the magnetic equator it is possible to consider 
the gravitational and nonlinear electrodynamic bending of 
rays as occurring approximately in the planes containing 
the neutron star center, the output aperture can also be 
presented as the sector of a ring of radius, p, width, dp., 
and angle opening 59. Thus the output aperture area will 
be equal to: S ou t = pdpS9. 
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Fig. 3. The HE dependence 8tp\ versus h/b for different values B 2 (b)/B 2 : (a) B 2 (b)/B 2 = 0.25, (b) B 2 (b)/B 2 
0.04, (c) B 2 (b)/B 2 = 0.01 (B 2 (6) is the square of magnetic field at the distance b from the neutron star). 



( 



f 



Fig. 4. The HE dependence 6<p* 2 versus h/b for different values B 2 (b)/B 2 : (a) B 2 (b)/B 2 = 0.25, (b) B 2 (b)/B 2 
0.04, (c) B 2 (b)/B 2 = 0.01 (B 2 (b) is the square of magnetic field at the distance b from the neutron star). 



Since the energy value transferred per unit time along 
the bunch of rays does not depend on the distance, then 
from the equality Ii n Si n — I utS 0U t we obtain: 



lout — 



h[ p 2 + {h+hf] bdb 
(If-R 2 ) pdp 



(17) 



In this relation values p and b are not independent, but are 
connected by the equation, which can easily be obtained 
from geometrical consideration (see Fig. lb): 

/ 2r \ 
p = bcosa + [6sina + l 2 ]tg[a - Sf* 12 j ■> ( 18 ) 

where ao = arcsin(&/li), r g is the Schwarzschild radius of 
neutron star, and 5ip\ 2 IS defined by espression (16). 

This equation is transcendental relatively to the im- 
pact parameter b. Its analytical approximate solution can 
be obtained only under a large number of restrictions on 
the parameters contained in it. In particular, if R ~ 10 2 
km, h - 10 4 km, r g - 0.2 km, l 2 - 10 17 km, Sip* 12 < 10~ 2 , 
then is transcendental equation (18) by 2R > b > R can 
be approximate by algebraic equation of the 7th order rel- 



atively to the impact parameter b : 



\b 2r a 



157r£77i i2 |m| ; 



h b 4b 6 

Since b> R, and R/l\ ~ 10 -2 , the last two terms in the 
squared brackets give the small correction to b/l\, which 
decreases with the increasing b. As a result, the approxi- 
mate solution of this equation can be written as: 



hP( 1+ h 
h L P 



15<77 l!2 |m| 2 Z 



2/6 



-0 



S l 2 



},(19) 



hp ilfp 6 \l\p 2 . 

and this solution is valid only by variation of p 
within the limits: 2l 2 R/h > P > h[R/h - 2r g /R - 
15<r7i, 2 |m| 2 /(4.R 6 )]. 

Substituting expression (19) in the equality (17), we 
then obtain, that in this field of variation of p the detected 
gamma-ray flux will be equal: 



lout — Iq 



1 - 



757r£?7 



1.2 



m\ 2 a ^r 2 J 3 



Alfp 7 



O 



iIp 3 



(20) 



Thus, gravitation and nonlinear electrodynamic bending 
of gamma-rays leads to a decrease of the detected flux in 
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the considered part of space. It is caused by consideration, 
that in the absence of the ray bending the illuminated area 
P > Rh/h is less than illuminated area p > ^[R/h — 
2r g /R — 157r^77 1)2 |m| 2 /(4i? 6 )] in the presence of bending. 
As a result of the ray bending, part of the gamma-ray 
energy flux is transferred from the field of space p > RI2/I1 
to the field of space bounded by the conic surfaces p = 
Rl 2 /h and p = l 2 [R/h - 2r g /R- 157r^i, 2 |m| 2 /(4 J R 6 )], 
thus decreasing the flux of energy outside the cone p > 
Rh/h- 

We will make some estimates. If we use the above pre- 
sented numerical values of R, h , r g , I2 and if to assume that 
in the case of some gamma-ray pulsars near the neutron 
star surface B = 2 • 10 13 G, (Thompson 2000) we obtain, 
that relative value (To — Iout)/Io changes from 10 -2 at 
p = Rl 2 /h ~ 10 15 km to 10~ 4 at p = 2Rl 2 /h ~ 2 • 10 15 
km. 

Due to the larger magnetic field Bo <~ -10 15 — 10 16 G, 
(Zhang & Harding 2000; Duncan & Thompson 1992; 
Thompson & Duncan 1995, 1996) in the case of a magne- 
tar this value will change in wider boundaries. However, 
the correct estimates in post-Maxwellian approximation 
can be made only in the frame of the BI electrodynam- 
ics, since in the HE theory, if the magnetic field is greater 
than the B q value, the Lagrangian should be changed dras- 
tically, because it sh ould contain logarithmic terms. 

5. Analysis of x-ray and gamma-ray astronomy 
technique applications for observing nonlinear 
electrodynamic effects 

The modern accuracy of gamma ray flux parameter mea- 
surements provided by the use of extra- Terrestrial tech- 
niques is much worse than the accuracy of measurements 
of similar parameters in the optical range (Boyarchuk et 
al. 1999). However, due to the opacity of the strongly mag- 
netized neutron star magnetosphere to the optical emis- 
sion, X-rays and gamma rays give the unique possibility to 
research for nonlinear electrodynamic effects in the strong 
magnetic field of a neutron star. 

5.1. The effect of gamma-ray flux dispersion 

Taking into account the capabilities of modern X-ray 
and gamma ray space observatories discussed in the 
Introduction, we will analyze, which of the above men- 
tioned effects can be observed using space astronomy tech- 
niques. 

In the case, when the distance l\ between the gamma 
ray source and the neutron star is comparable with the 
neutron star radius, the energy flux (20) is close to one 
even for the impact distances comparable to the star ra- 
dius. Thus, the rays which underwent significant nonlinear 
electromagnetic influence from the neutron star magnetic 
field can be weakly dispersed, and as a result their in- 
tensity at the point of observation will be sufficient for 
detection. 



Let us estimate the magnitudes of these effects. 
However, first of all it is necessary to make the following 
clarification. All the above discussed effects depend not 
only on the magnetic field value, but also on the choice of 
the model nonlinear vacuum electrodynamics (since tran- 
sition from one model to another changes the value of 
post-Maxwell parameters 771 and 772). At present it is im- 
possible to say, which value of these parameters is agrees 
with observational data. 

It is necessary to note, that in the HE nonlinear 
electrodynamics (Ritus 1986) the expansion parameter is 
p = (H/B q ) 2 . By B = 2 • 10 13 G this parameter is equal to 
p = 0.2 and post-Maxwellian expansion (5) of QED still 
can be used for estimation of nonlinear effects values. 

Thus, when making the estimates we will take into 
account the most well known models of nonlinear vac- 
uum electrodynamics, such as the BI electrodynamics and 
the nonlinear electrodynamics, which is the direct conse- 
quence of quantum electrodynamics. According to these 
theories the post-Maxwell corrections to the Lagrangian 
of nonlinear electrodynamics are about ~ 10~ 4 B 2 /£? 2 . 

However, it is necessary to note, that these estimates 
are quite conditional, and should be considered as esti- 
mates of concrete theories. It may be, that some other 
theory is more adequate to nature, and will give other 
estimates of the studied nonlinear electrodynamic effect 
values in the discussed experiments. During the space ex- 
periments it is necessary to search for and measure record 
all the above mentioned effects in order to obtain the val- 
ues of post-Maxwell parameters 771 and 772 from the results 
of observational data processing. Thus, the purpose of as- 
trophysical observations is not only the testing the predic- 
tions of one or another nonlinear vacuum electrodynamics 
model, but also what is more important the measuring of 
parameters rji and r/2- 

Using equalities (16), it is not difficult to obtain, that 
the maximum value of the nonlinear electrodynamic bend- 
ing angle of gamma-rays in the magnetic field of a pulsar 
with B = 2 • 10 13 G is about Sip* ~ 2.1 • 10~ 4 rad w 40 
arcsec according to the quantum electrodynamics predic- 
tions and Sip* ~ 1.4 • 10~ 5 rad w 2.8 arcsec according to 
the BI electrodynamics. 

In the field of gamma-ray pulsar or magnetar with B ~ 
10 16 G the maximum value of this angle increases: 5ip* ~ 
3.5 rad w 200° according to the BI electrodynamics. 

It is necessary to note that, the magnetic field of real 
neutron stars is most certainly not dipolar (Feroci et al. 
2001). In the case of a non dipolar field the effective 
volume of nonlinear electrodynamics interaction, will be 
smaller in comparison with the case of a dipolar field. 
Thus, the gamma-ray bending will be less clearly pro- 
nounced. However, the above obtained estimates made for 
the dipole field can be used for model quantitative evalu- 
ations of nonlinear electrodynamics effects. 

It is necessary to note that the obtained values char- 
acterize deflection angles at the source. At the point of 
observation possible deflection of a ray will be determined 
by the attitude Sb/h, were 5b is the variation of the im- 
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pact parameter of a ray because of the proper move of 
the source and I2 is the distance from the source to ob- 
server. For typical values l 2 = 10 17 km (3 kps), 8b = 10 2 
km, Sb/l 2 ~ 3.4 • 10 -12 arcsec that is beyond any mea- 
surable limits now. The only nonlinear electrodynamics 
effect, which can be measured principally, is the effect of 
gamma ray flux dispersion by the neutron star magnetic 
field. It follows from equation (20) that flux attenuation 
coefficient depends on the impact parameter b as well as on 
the ray bending angle. Thus, studying this effect we could 
also obtain information on the nonlinear electrodynamics 
bending of a ray in the source. 

The effect of gamma ray flux dispersion can be in- 
terpreted more accurately, if the neutron star and the 
gamma ray source are moving relatively each other with 
periodical time dependence of their mutual location. As 
our estimates show, the limits of the flux attenuation co- 
efficient value in this case are very wide. Hence, in the 
case of a rotating neutron star, if gamma rays are emitted 
in its vicinity, we will have periodic variations of the im- 
pact parameter which lead to the periodic variations of the 
flux attenuation coefficient. Such variations will efficiently 
modulate the mean light curve, which can be detected by 
the outside observer. Under certain conditions the same 
effect in the orbital light curves could be observed in the 
case of pulsars in binary systems. However, the presence of 
accreting matter can lead to some difficulties in resolving 
the pure nonlinear electrodynamic effect from the mean 
orbital light curves. 

5.2. The effect of birefringence 

The main qualitative difference in the predictions of 
different nonlinear electrodynamics theories, which can be 
observed in the considered astrophysical conditions, is the 
absence of birefringence in the BI theory. According to the 
HE theory vacuum birefringence effects the polarization of 
the emitted photons as well as the radiative opacity of the 
near neutron star medium. Such effects leads to a signa- 
ture in the emission spectra. However, for X-ray frequen- 
cies and pulsar magnetic fields to be quite pronounced 
this effect requires rather high plasma density (10 23 cm~ 3 ) 
(Meszaros & Ventura 1979), which can take place in some 
accretion powered pulsars. In the case of individual mag- 
netic neutron stars (gamma-ray pulsars and magnetars) 
the magnetosphere is quite transparent for hard emission. 
The main effect, which has influence on the propagation 
of emitted photons, is absorption by pair production in a 
strong magnetic field. This suppresses the radiation above 
about 1 MeV and effects the gamma ray beaming (Riffcrt 
et al. 1989). 

As for X- ray frequencies in the case of pulsars with 
poor magnetosphere we suppose, that the only observable 
difference between predictions of BI and HE theories is 
that according to the first theory the nonlinear electrody- 
namic bending angle in the dipole magnetic field of a neu- 
tron star and consequently the flux attenuation coefficient 



will be the same for any polarization of electromagnetic 
waves. In the HE theory r)\ 7^ rj 2 as well as the nonlin- 
ear electromagnetic bending angle and consequently the 
flux attenuation coefficient depend on the electromagnetic 
wave polarization. 

Thus, this theory predicts that if the gamma-ray 
source moves periodically relatively to the neutron star 
(i.e. in the case of rotating neutron star or pulsar in bi- 
nary system) we will obtain different mean light curves for 
emitted electromagnetic waves with perpendicular mutual 
linear polarization. 
6. Discussion 

Although the observation of the manifestations of the non- 
linear electrodynamics effects in astrophysical objects re- 
quires special conditions, in principle, they can be ob- 
served. The main astrophysical objects, where the non- 
linear electrodynamic effects can be revealed more clearly, 
are certain kinds of gamma-ray pulsars and magnetars. 
These effects can be manifested as some peculiarities in 
the form of their hard emission pulsation. 

The typical luminosity of magnetars and certain kinds 
of rotation-powered pulsars in hard emission is about 
10 34 _ 1Q 36 erg s -i (Mereghetti 2000). For example, the 
sensitivity level of the Chandra instruments for 10 4 s ex- 
posure corresponds to the luminosity of galactic objects 
- 10 30 - 10 31 erg s^ 1 ( Garcia et al. 2001). Hence, for 
the exposure time of about 1 s Chandra permits to obtain 
the mean light curve of a source with the luminosity of 
<~ 10 32 — 10 33 erg s -1 , this is about two order lower than 
the typical luminosity of such magnetar-candidate sources 
as AXPs or SGRs. 

The imaging mode of an instruments is not necessary 
for obtaining the mean light curve, because for most of 
these objects (soft gamma-ray repeaters, anomalous X- 
ray pulsars, etc) the pulsation period is known. However, 
the detector outputs should be folded over the maximum 
possible time of continuous exposure with that known pe- 
riod. For example, we can take 10 5 s as the exposure time 
of a given source, thus according to the announced sensi- 
tivity for TeCd detector of the INTEGRAL IBIS instru- 
ment (Winkler 2001) we obtain the minimum detectable 
intensity of periodic processes (at 8a level), which is about 
1 mCrab. This estimate permits to obtain the mean curve 
of pulsation for a source like a SGR even in its quiescent 
state. Since the number of known magnctar-likc sources 
is not large, the inspection of most of them during the 
mission will be quite justified. 

Measurements of the bending angle in the source need 
observations of at least the proper orbital motion of the 
pulsar, which is not accessible by currently operating in- 
struments. 

The additional group of effects may be connected with 
such objects as accretion-powered binary X-ray pulsars 
and, possibly other kinds of tight binaries containing mag- 
netic neutron stars. As it was mentioned above, the objects 
with most suitable conditions for nonlinear electrodynam- 
ics effects are binary systems with underfilled Roche lobe 
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supergiants as an optical companion. The value of the or- 
bital period provides the necessary time for continuous 
observations of these sources. Because most of pulsar sys- 
tems with underfilled Roche lobe optical companion are 
characterized by large orbital periods (4U1538-522, 3.73 
d; 4U1907+097, 8.38 d; 1E1145. 1-614, 5.648 d; Vela X- 
1, 8.965 d) (Bieldsten et al. 1997), rather long exposures 
are necessary to obtain detailed mean orbital light curves. 
The existense of certain non-typical forms can indicate 
the presence of nonlinear electrodynamic effects in such 
objects. 

Acknowledgements. We are grateful to Academician Georgiy 
T. Zatsepin for fruitful stimulating discussions. Authors would 
also like to express their warm thanks to Ekaterina D. Tolstaya 
for grammar corrections and Alexander A. Zubrilo, Vitaly V. 
Bogomolov, Oleg V. Morozov for the help with preparing this 
manuscript. 

Part of this work was supported by the Russian Foundation 
of Basic Research. 

Appendix A: Eikonal equation 

Using Lagrangian (9), it is not difficult to obtain the ex- 
pansion of vectors D and H onto degrees £ with accuracy 
up to £ 2 inclusive: 

D = E + £{2??i(E 2 - B 2 )E + 4r/ 2 (B • E)B} + (A.l) 

+e 2 {3ry 3 (E 2 - B 2 ) 2 E + m (B • E) 2 E+ 

+r ?4 (E 2 -B 2 )(B-E)B}}, 

H = B + £{2r/i(E 2 - B 2 )B - 4r/ 2 (B • E)E}+ 

+£ 2 {377 3 (E 2 - B 2 ) 2 B + m (B • E) 2 B- 

-7? 4 (E 2 -B 2 )(B-E)E}. 

Under the assumption of a "weak" plane electromagnetic 
wave (8) with the use of an approximation linear in vectors 
b and e, we can obtain from relation (Al) and equations 
(1) a uniform system of three linear algebraic equations 
relatively to three e components: 

II a % = 0, (A.2) 
where 

U] 2 

n af} = A Q {k a k + k 2 )S a0 }+ 

and for more compact record we introduce the following 
notations: 

Aq = 1- 2 ??1 £B 2 + 3 Vs eK M = - 12%? 2 Bg, 

A 2 = 4r? 2 £ - mt 2 K , N B = [k x B ] . 

For the existence of nontrivial solutions of the equations 
system (A2) it is necessary, that det\ |n Q/3 | | = 0. The deter- 
minant of tensor H a P can be calculated in the most simple 
way, if we use the formulas of tensor algebra verified in the 
work of Denisova & Mehta (1997). 



The condition of equality to zero of the second order 
tensor determinant in three-dimensional Euclidean space 
H al3 can be rewritten using these formulas in the form: 

2n (3) -3n (1) n (2) + nf 1) =o, 

where 

n(jv) = n aiJ g 1 (5 l 3 ia2 ii a2J 3 2 (j / 3 2a 3 . . . u aN p N 5p NCtl . 

To compose the tensor's degrees n Q/3 , after the reduction 
by uj 2 A /c 2 we obtain the following dispersion equation: 

^{a 2 A B 2 + Al] + ^{A 4 [k x B ] 2 B 2 - (A.3) 
-A 2 A B 2 k 2 - 2A 2 k 2 + ArA [k x B ] 2 - 
-A 2 A (k x B ) 2 } + A 2 X - AiA k 2 [k x B ] 2 + 

+A 2 A k 2 (k • B ) 2 - A 4 [k x B ] 2 (k • B ) 2 = 0, 

where A 4 = A X A 2 - A\. 

We will find the solution lu — uj(k) of this equation as 
the expansion onto degrees £, similar to the expansion of 
the Lagrangian (9): 

LU = ck[l + iF + i 2 U], (A.4) 

where F and U are unknown functions. 

Let us substitute now this expansion into equation 
(A3). Because the Lagrangian (9) is written with accu- 
racy up to terms proportional £ 2 inclusive, after calcula- 
tions we can neglect the insignificant terms £ 3 . As a result 
we obtain a relation, which has the form of expansion in 
degrees £, in view of its complication we will not write it 
here. 

To make equal to zero the coefficients of this expansion 
in degrees £, we obtain the equation allowing to determine 
the unknown functions F and U. In the lowest approxima- 
tion we have: 

{Fk 2 + 277! [k x B ] 2 }{Ffc 2 + 2 V2 [k x B ] 2 } = 0. 

To resolve this equations relative to F and substitute it 
by (A4), and obtain the dispersion equations (10). 

Appendix B: Calculation of the ray bending angle 

Differentiating the eikonal (15) with respect to a and mak- 
ing it equal to the constant /?, we obtain equation for a 

ray ip = <p(r) : 

r 

/dr 
; - (B.l) 
r 2 y/§[l + ^]-# 

If we denote the source emission frequency as ujq, then 
the constant £ = w , and a = —uj b/c, where b - is the 
impact distance of the ray. 

As a result expression (Bl) becomes: 

r 

/dr 
/ (B.2) 
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The integral in the right-hand part of this equality can be 
expressed through elliptical functions. However, it is more 
favorable for our purposes to find another approach. Let 
us differentiate the equality (B2) with respect to tp and 
then to raise it into the second-order power: 

4?7i£m 2 6 2 " 



{ dry 



r 



(B.3) 



To find the solution of this equation we will use the well 
known Darwin method (1961). First of all, we will intro- 
duce the subsidiary variable u = 1/r. The equation (B3) 
becomes: 



(|) 2 = p[ 1+4 ^ mV - tV l- 



(B.4) 



We will find the solution of this equation in the form, as 



u = A + Qi sin*(<^), 



(B.5) 



where A and Qi are constants. 

Substituting expression (B5) into equation (B4) and 
restricting ourselves by the accuracy linear in the small 
parameter ?7i£m 2 M 6 , we obtain: 

A.O, 0l = i[l + ^], 

and function ^{tp) with acceptable accuracy should satisfy 
equation 

d*(p) , 2r/ 1 £m 2 



= 1 - 



dip 6 6 
The solution of this equation is: 
77i £m 2 



[3 — 3 cos 2 tp + cos 4 tp] 



*(» = ip - ipa 



166 6 



[60(p — 16 sin 2(p + sin &<p\ 



In the case, when the gamma-ray source is located at a 
limited distance from a neutron star or even in its nearest 
vicinity let us consider that this source is at the point 
r = Zi, tp = it. Then for a ray with impact distance b the 
solution of equation (B4) can be represented as: 

77i£m 2 



u = — = Qi sin {ip 

r 



-(po- 



rn 6 



■ [60</?— 16 sin 2(/3+sin &tp\ } . 



Because this ray should pass through the point r = h, tp = 
it, the integration constant tpo is more complicated in com- 
parison with expression (B6): 



157rr7i^m 2 
466 



• / 1 X 

+ arcsin(— -) 
Win 



The bending angle 8tp\ of a ray after its passing through 
the neutron star magnetic field will be equal to: 5tp* = 
<P>i — Pd, where tpi is the angle of a ray inclination to the 
X axis at the point r = h, tp = w, and tp D is the angle of 
the detected ray inclination to the X axis. 

Since in the considered case I2 » h, the tp D angle 
can be determined from the condition: u(tp D ) = ~ 0. 
It follows from this, that 



•PD 



157r?7i^m 2 
466 



+ 



1 + 



15?7i^m ; 



2-, 



46 6 



1 



arcsin( )+ 



Q 



r/i^m 2 r 1 1 , 

+ ^666-^ Sin[4arCSin( Q^ )] " 16sin[2arcsin(^-)]}. 

As it is well known, the tangens of tpi inclination angle at 
the point r = h, tp = w, is equal to the derivative dy/dx 
at this point: 

dy , 



tg(<Pi) = -j-{r = h,tp = 7r) 



r' sin tp + r cos tp 



1 



7 Op = 7r )' 



r' cos tp — r sin tp l\u' 

where the apostrophe denotes the derivative over tp. 

In the first approximation in the small parameter 
£m 2 /6 6 it can be obtained from this: 



tpi = arctan 



bQi 



in 



As a result we obtain the formula (16) for the bending 
angle 5tp\ of a ray. 
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